Astronomy & Astrophysics manuscript no. rxj 1347 


© ESO 2008 


February 2, 2008 





The mass distribution of RX J1 347-1 145 from strong lensing'*' 

A. Halkola^**, H. Hildebrandt^ T. Schrabback^ M. Lombardi^, M. Bradac^'^, 
T. Erben\ P. Schneider^ and D. Wuttke^ 

^ Argelander-Institut fiir Astronomic, Auf dem Hiigel 71, 53121 Bonn, Germany 

^ European Southern Observatory, Karl-Schwarzschild-Strasse 2 85748 Garching bei Munchen, Germany 

^ Kavli Institute for Particle Astrophysics and Cosmology,2575 Sand Hill Rd. MS29, Menlo Park, CA 94025, USA 

^ Department of Physics, University of California, Santa Barbara, CA 93106, USA 

Received October 2007 / Accepted January 2008 

ABSTRACT 

Aims. We determine the central mass distribution of galaxy cluster RX Jl 347- 1145 using strong gravitational lensing. 
Methods. High resolution HST/ACS images of the galaxy cluster RX Jl 347- 1145 have enabled us to identify several new multiple 
image candidates in the cluster, including a 5 image system with a central image. The multiple images allow us to construct an accu- 
rate 2-dimensional mass map of the central part of the cluster. The modelling of the cluster mass includes the most prominent cluster 
galaxies modelled as truncated isothermal spheres and a smooth halo component that is described with 2 parametric profiles. The 
mass reconstruction is done using a Markov chain Monte Carlo method that provides us with a total projected mass density as well as 
estimates for the parameters of interest and their respective errors. 

Results. Inside the Einstein radius of the cluster (~35'', or ~200 kpc, for a source at redshift 1.8) we obtain a total mass of 
(2.6 + 0.1) X 10^"^ Mq. The mass profile of the cluster is well fitted by both a Navarro, Frenk and White profile with a moderate 
concentration of c = 5.3^^^ and r2oo = 3.3^^^ Mpc, or a non- singular isothermal sphere with velocity dispersion cr = 1949 + 40 km/s 
and a core radius of r^. = 20 + 2 The mass profile is in reasonable agreement with previous mass estimates based on the X-ray emis- 
sion from the hot intra-cluster gas, however the X-ray mass estimates are systematically lower than what we obtain with gravitational 
lensing. 
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1. Introduction 

A precise knowledge of the masses and mass profiles of galaxy 
clusters is a key to better understand what the mass and energy in 
the universe is made of. The standard model of cosmology that 
■ has emerged over the past 100 years is remarkably successful 
I in reproducing a wealth of observational phenomena. These in- 
, elude the formation and evolution of structure from the tiny tem- 
perature variations seen in the cosmic microwave background to 
^ galaxies and clusters of galaxies, the element abundances pro- 
, duced in the Big Bang nucleosynthesis and the apparent accel- 
eration of the expansion of the universe as seen from the bright- 
nesses of distant supernovae. The greatest drawback of the stan- 
dard model is the need to include two dark components in the 
energy density of the Universe. The first is dark matter which is 
needed to explain the formation of structure and the dynamics of 
both galaxies and clusters of galaxies. The second is dark energy 
which in turn is needed to make the Universe flat and to explain 
the supernova data. Both of these dark components are still un- 
explained although candidate particles for dark matter exist. 

Clusters of galaxies are of prime interest when trying to con- 
strain the properties of both dark matter and dark energy. The 
mass budget and therefore the dynamics of clusters are dom- 
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inated by dark matter, whereas the mass function and power 
spectrum of clusters is influenced by both dark matter and dark 
energy. In order to obtain strong constraints on the various cos- 
mological parameters such as Dm, os and the equation of state 
w of dark energy, accurate mass estimates for many clusters are 
necessary. Gravitational lensing is a powerful tool for determin- 
ing cluster masses since the deflection of light is independent 
of the nature and dynamical state of the matter in clusters and 
can therefore provide unbiased estimates of the total masses of 
galaxy clusters. 

An important prediction of numerical simulations with cold 
dark matter is the so called universal mass profile. For large radii 
(greater than the scale radius) the density p depends on r like 
p(r) oc r"^, the exact value of the slope at small radii (inside 
the scale radius) is unknown but is believed to be somewhere 
between -1.5 and -1. Although there is strong support for the 
universal mass profile, other mass profiles such as the family of 
isothermal sphere profiles have not been ruled out. By studying 
the mass profiles of galaxy clusters we can potentially constrain 
the nature of dark matter through comparison with simulations. 

As the most luminous X-ray cluster (ISchindler et al.lll995l) 
known to date RX J1347-1145 is obviously of great inter- 
est. It has been studied spectroscopically by ICohen & KneibI 
(i200Z) and exten sively with X-rays and Sunyaev-Zel'dovich 
eflTect (Hc hindler efaD 119971: IPo intecouteau et al. 19991 I20Q1I: 
lAflen etal.1 i2002t iGitti & Schindlen »2004c iGitti et alTl20Q7l) . 
Gravitational lensing based only on the arcs has also been used 



to estimate the mass ^Sahu et al. 1998) as well as weak lens 



ing (iFischer & Ty sonlfl 997:. Klin g et al.ll2005l) . Recently a com- 
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bination of both the weak and strong lensing has also been used 
(iBradac et al.ll2Q05l) . The lack of space-based observations with 
the Hubble Space Telescope (HST) has so far limited the number 
of multiple images available for the strong lensing modelling. In 
this paper we use deep HST images taken with the Advanced 
Camera for Surveys (ACS) to identify new multiple image candi- 
dates that enable us to derive well constrained mass maps for the 
central regions of RX Jl 347- 1145 . We have obtained spectra 
of some of the multiple images with the F0RS2 spectrograph on 
the Very Large Telescope in Chile. The spectroscopic redshifts 
obtained are crucial in fixing the mass scale of the cluster. 

The cosmology used throughout this paper is Dm =0-30, 
Qa=0.70 and Ho=70 km/s/Mpc, unless otherwise stated. With 
this cosmology an arcsecond at the redshift of the cluster 
(Zci=0.451) corresponds to 5.77 kpc. 

2. Data 

2.1. HST Imaging 

Strong gravitational lensing requires deep, high quality and high 
resolution images to facilitate the identification of multiple im- 
ages in galaxy clusters. The best instrument for this is the 
ACS/WFC on bo ard the HST. Follow ing the strong and weak 
lensing analysis of'Brada c et al.l (l2005l) . we have acquired imag- 
ing data in filter bands F475W, F814W and F850LP, 5280s in 
each band (PI: T. Erben, proposal no. 10492). This gives us not 
only high resolution images but also limited colour information 
for the galaxies. 

The data reduction is based on the bias and flat-field 
corrected images obtained with the ACS calibration pipeline 
CALACS. Before the standard pipeline reduction proceeds we 
compute a noise model for each exposure. These are later used 
to appropriately weight the individual exposures when coadding 
the exposures. At this point also the badpixel masks are updated. 
For more details on the noise modelling and bad pixel masking 
see Marshall et al. (in preparation). 

The two CCDs on the ACS have 2 readout ports each, and 
this can sometimes lead to residual bias level between the four 
image quadrants. We have therefore subtracted the sky back- 
ground separately in each of the fou r image quadrants . 

We finally use MultiDrizzle (iKoekemoer et al.1 l2002h to 
coadd all the individual exposures of each filter. This step also 
corrects for geometric distortion caused by the telescope optics 
and detectors, and rejects the cosmic rays present in individual 
frames. In order to accurately align all the images taken in the 
3 bands, we select highS/N objects in the individual exposures 
and compare their positions to find out any residual shifts and 
rotations between the images. The necessary corrections given 
to MultiDrizzle are calculated using the geomap task in IRAK 
In the final coaddition we use the noise models calculated ear- 
lier to weight the individual exposures using inverse variance 
weighting. We apply the "square" drizzle kernel in MultiDrizzle 
in combination with a reduced pixel scale of 0^.^03 and slightly 
shrunk pixels (pixfrac=0.9). More details on the data reduction 
can be found in Marshall et al. (in prepararion). 

2.2. Photometric redshifts 

The photometric redshifts used in this work are calculated from 
the 5 band catalogue based on images taken with the VLT/FORS 
in the U, B, V, R, an d / bands. For mo re details on the data 
and data reduction see lBradac et al.l (|2005|) .The photo metric red- 
shifts are calculated using the Hyperz package (.Bolzonella et al.l 



l2QQQh . We do not use any prior on the magnitude which sup- 
presses the probability density at high redshifts. The multiple 
images we are interested in are often highly magnified which 
increases the flux coming from a background galaxy. In many 
cases the multiple images are still very faint however and it is 
very difl&cult to obtain good photometric redshifts for the im- 
ages. The three ACS bands on the other hand are not enough 
to give a photometric redshift estimate on their own, and com- 
bining them with the ground based date can introduce system- 
atics/biases. To check this, we have also includes the ISAAC 
^^-band data in the photometric redshift estimation. The best-fit 
photo-z's of the multiple images are barely aff'ected and addition- 
ally some of the probability distributions even get broader and 
less well defined when using the ^^-band. This might be due to 
slight photometric miscalibrations and difficulties in the equali- 
sation of the point spread function of the diff'erent images leading 
to systematic colour errors. The photometric redshifts probabil- 
ity distributions are compared to the lensing ones in Sects. 15.31 
and [5A2l 

2.3. VLT Spectroscopy 

Spectra were taken for potential cluster members, a magnitude 
limited sample of galaxies (to be used as a calibration sample 
for photometric redshifts) as well as multiple image candidates. 
The spectra were taken with the F0RS2 spectrograph at the 
Very Large Telescope of the European Southern Observatory. 
Reduction of the spectra and detailed results will be presented 
in a forthcoming publication (Lombardi et al. 2008, in prepa- 
ration). In this work we only use the spectra and spectroscopic 
redshifts obtained for multiple image systems as detailed in later 
sections. 



3. Components of the Strong Lensing Models 

In this section we describe the strong lensing modelling of the 
cluster in detail. 

3.1. Mass components 

The masses of clusters of galaxies have three major (distinct) 
components: the galaxies and their DM haloes, the hot intra clus- 
ter gas and the dark matter (DM) halo of the cluster. Although 
the galaxies contain only a small fraction of the total mass they 
have a significant local eff'ect on the multiple images and need 
to be included in the modelling. The galaxies in a cluster aff'ect 
strongly the formation of strong gravitational lensing features 
such as arcs and multiple images (e.g. Meneghetti et al. 20Q7|). 
It has also been investigated how the increased cross section of 
galaxy- sized haloes in denser group or cluster environments can 
potentially be used to constrain the mass profiles of cluste rs and 
the dark matter content of the galaxies themselves (lKingll20()7l: 
Tu et al. 2007). In this work, however, we do not attempt to do 
this. We assume that the hot X-ray gas and the DM together can 
be described by two dark matter haloes. 

3.1 .1 . Cluster Galaxies 

We have identified the cluster galaxies from the cluster red se- 
quence. The red sequence in the ACS F475W-F850LP, F814W 
colour-magnitude diagram is shown in Fig.[T] We have selected 
lensing galaxies based on their magnitude and F475W-F850LP 
colour. All galaxies brighter than 23 (AB) in the F814W pass- 
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Fig. 1. The red sequence of galaxy clusters in the 
F475W-F850LP, F814W colour-magnitude diagram. The 
galaxies that are included in the mass model are marked 
with red circles and have F475W-F850LP colour within 0.5 
magnitudes from the mean of the red sequence. 



band with F475W-F850LP colour within 0.5 magnitudes from 
the mean of the red sequence are included. We have additionally 
included all spe c trosc opically confirmed cluster galaxies from 
ICohen & Kneibl (12002). This results in a total of 119 galaxy 
lenses in our lensing modelling. 

The profiles of the cluster gala xies are based on a truncated 
isothermal sphere profile ( Brainerd et al.| [T996). The profile is 
characterised by a velocity dispersion cr and a truncation ra- 
dius s. The 3 dimensional mass density of a truncated isothermal 
sphere is given by: 



Phhs(r) = 



(1) 



The galaxy haloes can be elliptical with the ellipticity im- 
pleme nted following the method presented in (iGolse & Kneibl 
I2002h . All the galaxies included are treated in exactly the same 
way. This includes also the two massive central galaxies and the 
galaxies near multiple images. 

The velocity dispersions of the galaxies are estimated from 
their F814W magnitude using the Faber- Jackson relation. Since 
we have been unable to find a good reference point to a fun- 
damental plane or a Faber- Jackson relation at the redshift of the 
cluster we have used the fun damental plan e kno wn for the galaxy 
clusters A2218 and A2390 (IZiegler et al. 11200 iHFritz et al.ll2005l) 
at redshifts 0.18 and 0.23 respectively. These same clusters were 
used as reference points when estimating the velocity disper- 
sions of the galaxies in A 1689 in H alkola et al.l (I2006D . To model 
the galaxies in RX J 1347 we have taken a model spectrum of 
an early ty pe elliptical g alaxy with a formation redshift of 5 
from Bruz ual & CharlotI (12003 ^), and passively evolved it from 
the redshift of RX J1347 (z = 0.451) to 0.2. This model is used 
to calculate both the K-correction from the obser ved F814W to 
the F775W of the Faber- Jackson relation shown in lHalkola et al.l 



(l2006h for clusters at redshift 0.2, and the evolution of the galaxy 
colours. The measured F814W magnit udes a re corrected for 
galactic extinction according to Schlegell^lalJ ^99 8) obtained 
using the NASA/IPAC Extragalactic Database (NED). 

We assume a linear relation between the truncation radius s 
of a galaxy and its velocity dispersion cr. This is expected on the- 
oretical grounds (Merritt 1983). Recently this topic has also been 
investigate d using simulations that t ake into account the eff'ect of 
baryons by iLimousin et al.l (l2007bl) . These N-body simulations 
were done in order to study the tidal stripping of dark matter 
haloes of galaxies in a cluster, including star formation so that 
the stellar components of the galaxies can be followed. The sim- 
ulations agree with tren ds observed using galax y-galaxy lensing 



ulations agree with tren ds observed using galax y-g; 
in galaxy clusters, e.g. iNataraian et al] (2002h: Li 



Limousin et al.l 



(2007a). It was additionally shown in lHalkola et al.i (i2007i) that 
for the strong lensing analysis of Abell 1689 (which has much 
stronger constraints on the mass from the multiple image sys- 
tems identified in the cluster) it was not possible to constrain 
the slope of the power-law relation between s and cr. Since the 
multiple images in RX J 1347 do not constrain the sizes of the 
galaxies significantly we assume the same truncation law for the 
galaxies in R X J1347 as was found for galaxies in the galaxy 
cluster A1689 (iHalkola et al.|[2Qb7i) as the clusters are both sim- 
ilarly massive. 

3.1 .2. Smooth cluster mass 

The X-ray measurements as well as the combined strong and 
weak lensing analysis of the cluster indicate that there is possi- 
bly a mass extension to the south-east of the cluster centre. The 
mass maps are otherwise fairly smooth and we therefore model 
the smooth mass component of the cluster with two paramet- 
ric profiles. The second halo is also necessary in order to accu- 
rately reproduce the multiple images observed. Using paramet- 
ric haloes to reconstruct the mass of galaxy clusters, as opposed 
to pixelated mass reconstructions for example, have several ad- 
vantages but also disadvantages. The simple description of the 
mass distribution gives generally very good convergence prop- 
erties when optimising the model parameters. On the other hand 
the mass profile is also relatively restricted in the mass distribu- 
tions that can be recovered. Although the second halo could in 
principle be connected to a physical structure in the cluster, it is 
also possible that a secondary halo is necessary in the modelling 
in order to accurately model the mass distribution of the cluster 
with only one distinct halo. These two haloes are both des cribed 
by either Navarro-Frenk- White (NFW Navarr o et aO 19971) or by 
non-singular isothermal ellipsoid (NSIE) haloes. 

The NFW halo is a prediction of cosmological numerical 
simulations and its 3 -dimensional mass density is given by 



Pnfw(^) = 



Pc 



(r/rs)(l + r/rs)2' 



(2) 



where pc is the characteristic density of the cluster and is the 
scale radius which marks the transition from p oc r"^ to p oc r"^. 

The NSIE profile is based on the popular isothermal sphere 
profile with the inclusion of a constant density core. The 3 di- 
mensional mass density of the NSIE profile is given by 



Pnsie(^) = 



1 



InG (r^ + rl) ' 



(3) 



where cr is the velocity dispersion of the isothermal sphere and 
Tc is the core radius of the sphere. 
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Both mass profiles have 6 free parameters: their position on 
the sky (a, S), the ratio of their semi-minor and semi-major axes 
r = b/a, the position angle 6, and the two parameters described 
above that define the density profile. 

3.2. Multiple images 

In the literature there are several published multiple imag e can- 
didates in RX J1347 fe.g. lSahu et aDll998l:lBradac et al.l l2QQ5). 
In addition to these we have identified several new ones in the 
HST ACS images of the cluster, bringing the total number of 
candidate lensing system to 13. The images are shown in Fig.[2l 
and some of their properties are shown in Tab. lA.ll and thumb- 
nail images in Figs. lA.ll - IA.6I of the appendix. The more ten- 
tative image candidates in a system are marked with a question 
mark in Tab. lA. II and Fig.[2l The following sections discuss the 
multiple image systems in some detail. 

The main analysis in this paper is based on image systems 1 
and 2 only. The locations of the 5 images in image system 1 pro- 
vide strong constraints on the distribution of mass in the cluster 
while the spectroscopic redshift of image system 2 is important 
in fixing the total mass scale. In section [54l we repeat the analy- 
sis with all image systems to check the eff'ect of including all of 
them, and to obtain lensing redshift probability distributions for 
the image systems. 



3.2.1. Image system 1 

This is a multiple image system with 5 images on 4 sides of the 
cluster and a central image. The identification is based on the 
similar colours and morphology of the images, all images have 
a brighter knot at one end. Also the photometric redshift proba- 
bility distributions of the multiple images in this system are very 
similar. They are flat and fairly broad and place the multiple im- 
age system in a redshift range 0.5-2.0. The Mf475w - Mfsmw 
and Mf814w - Mfssolp colours of the multiple images are also 
in good agreement. We have obtained spectra for all the multi- 
ple images, but have not been able to identify any emission lines 
in the individual spectra. The many skylines in the spectra also 
hinder us from cross correlating the individual spectra to con- 
firm the compatibility of the spectra of the multiple images. This 
mean that we cannot confirm the multiple image nature of the 
source from the spectra. Although the redshift cannot be fixed, 
the geometry of the images provide significant constraints on the 
mass profile. 



this galaxy is a background galaxy and must be lensed. In the ab- 
sence of counter images it cannot be eff'ectively used to constrain 
the mass profile. Also the lack of multiple images can constrain 
the mass distribution, see I Julio et al.l (l2007l) . This is only eff'ec- 
tive if the singly lensed images are relatively close to critical re- 
gions so that multiple images are formed for model parameters 
that fulfill the other constraints reasonably well. In our analysis 
we have only afterwards checked that no additional images for 
this system are predicted. 

3.2.4. Image system 4 

Similarly to image system 3 this one consists of only one dis- 
torted image and therefore cannot eff'ectively be used to con- 
strain the mass model. The spectroscopic redshift of this galaxy 
is Zspec = 0.785 (Ravindranath & Ho .200Z) . No additional im- 
ages for this system are predicted. 



3.2.5. Image system 5 

This image system consists of 3 faint images that are split by a 
cluster galaxy. The images are very sensitive to the mass in the 
galaxy which can therefore influence the total mass distribution 
significantly when this image system is included in the analysis. 



3.2.6. Image system 6 

A faint thin arc at a cluster centric radius of -12" larger than 
image system 2, and hence expected to be at a highe r redshift. No 
counter images found. This is image system C in iBradac et al.l 
(.i2005h . 



3.2.7. Image system 7 

An arc that nicely bends around the galaxy at Qf=13:47:34.8, 
^=-11:45:01.5. The exact shape of the arc is sensitive to the 
mass of the galaxy. 



3.2.8. Image system 8 

This is a 2-image syst em me r ging o n a tangential critical line, 
also identified in B radac et al.l (120051) as a possible arc candidate. 
There is a tentative third image on the east side of the cluster but 
a positive identification is difl&cult due to the faintness of the 
source. 



3.2.2. Image system 2 

Our image system 2 is composed of two ima ges. The brighter of 
the two arcs was already identified by Schindler et al.l (Il995h as 
the sou th-west arc. The fainter arc was later found by lSahu et al] 
('1998'). We have taken spectra of the two arcs and obtained a 
redshift of Zspec = 1.75 for the two of them; for more details see 
Lombardi et al. (2008, in prep.). The redshifts of the images fixes 
the overall mass scale of the cluster a nd provides s t rong s upport 
for them being multiple images. In iBradac et"aD (l200^ these 
correspond to images A4 and A5. Our models do not predict any 
additional detectable images for this system. 



3.2.9. Image system 9 

2 images merging on a radial critical line. We have identified a 
possible third image in the south-east part of the cluster. 

3.2.10. Image system 10 

A similar merging of 2 images at a radial critical line, close to 
image system 9. This one is much fainter though and less promi- 
nent and consequently less certain. A third image is expected 
near the third image of image system 9 but as this source is sig- 
nificantly fainter we have not identified one. 



3.2.3. Image system 3 

Image system 3 consists of only one di storted image. The s pec- 
troscopic redshift is Zspec = 0.806 (.Ravindranath & HqI|2002I) , so 



3.2.11. Image systems 11 and 12 

In the north-west of the cluster there are several images that de- 
fine a tangential critical line. The colours and morphology of the 
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Fig. 2. The multiple image systems identified in RX Jl 347- 1145 overlaid on a colour image of the cluster made from the images 
taken with the ACS in filter F475W, F814W, and F850LP. The box size is -880 kpc, north is up, east is left. Uncertain image 
candidates are marked with a question mark. Only image systems 1 and 2 considerably constrain the mass profile of the cluster. 



images vary, suggesting that the images originate from several 
diff'erent sources at a similar redshift. Possible counter images on 
the south-east and south west sides of the cluster are identified 
for system 11. The fainter image system 12 is more challenging 
in this respect althoug h several candidates have been identified. 
In lBradac et aP (l20Q5l) identified as D1-D4. 

3.2.12. Image system 13 

A bluish galaxy oriented tangentially with respect to the cluster 
centre. It is potentially lensed although we have not been able to 
identify any counter images to this one. 



4. Estimating the cluster mass distribution and its 
error 



In order to estimate the mass distribution of the cluster we use 
Bayesian statistics with Markov chain Monte Carlo methods to 
estimate the probability distribution function of the free param- 
eters of the cluster mass model. 

Our mass models have 13 free parameters in total, the red- 
shift of multiple image system 1 and 2x6 parameters for the 2 
smooth haloes for the global mass distribution. 
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4.0.13. Model;^^ 



Both the positional information as well as the relative magnifi- 
cations of the images in an image system are used in calculating 
the;^^ for the models. The;^^ is given by 



X 



N m 

2 



N Hi 

+ yy 



( - fujli^ujf 



(4) 

In equationlH N is the number of multiple image systems, ui the 
number of images in system / and fiij denotes the source posi- 
tion of image j in image system /, j^ij is the model magnification 
at the image position, fij and eij represent the isophotal flux and 
its estimated error of image j in image system /. The first dou- 
ble summation in equation [4] therefore represents the scatter of 
the source positions weighted by the model magnification. The 
errors are assumed to scale with the square root of the magni- 
fication since the magnification is given by the ratio of an area 
of the image in the image plane to that of the unlensed image in 
the source plane. It was shown in Hal kola et al.. (2006) that this 
is a good estimate for the;^^ calculated in the image plane. The 
second double summation includes the relative magnifications 
between the images of an image system in the modelling. This is 
calculated by comparing the unlensed source fluxes of the mul- 
tiple images. We use the first multiple image in each system as a 
flux reference. 



4.0.14. Bayesian statistics 

B ayes' theorem deals with conditional probability distributions. 
We are specifically interested in constraining a model with pa- 
rameters w given data d (and possibly some other prior knowl- 
edge of the system). The probability of a parameter set w given 
the data d can be written as 



P(w\d) ■ 



P(d\w)P(w) 
P(d) ' 



(5) 



where P(d\w) is the probability of the data d given the param- 
eters w, P(w) is the prior on the parameters w and P(d) is the 
probability of obtaining the data in the first place, also known 
as evidence. P(d\w) is measured by the of our models and 
we will assume that it is proportional to exp(-;^^/2). We as- 
sume flat priors P(w) for the parameters that are reasonable for 
DM haloes of galaxy clusters. P(d) is treated as a normalising 
constant. With these simplifying assumptions it only remains to 
estimate P(w\d). We do this with a Markov chain Monte Carlo 
method described below. 



2 1n(/7) = ^ 



, forx^(Wi+i) <X^(m) 



. (6) 



The new perturbed parameter sets w/+i are generated by random 
selection from Gaussian distributions centred on the individual 
parameters in w/, and the widths of the Gaussian distributions 
are matched so that when each parameter is varied individually 
the efl'ect on the acceptance rates converge to the same value of 
-60%. In the end all the widths are scaled so that an acceptance 
rate of -60% is achieved also when all the parameters are varied 
simultaneously. Another possibility would be to assume that the 
width is given by the error in the parameters. This would how- 
ever require prior knowledge on the errors of the parameters. 
By tuning the widths according to the acceptance rate we ensure 
that all free parameters of the model have the same weight in the 
Markov chain. The acceptance rate of 60% is chosen as a com- 
promise between denser local sampling of the parameter space 
at the cost of poorer global sampling (smaller width, higher ac- 
ceptance rate), and an attempt to sample quickly a large volume 
of the parameter space that is hindered by low acceptance. 

Each Markov chain is started with a mass model that has 
mass only in the galaxies. The two haloes describing the smooth 
mass component are set to be circular and to have no mass 
(cr = km/s for the NSIE profiles, and Pc = Mo/kpc^ for the 
ENFW profiles). The two haloes are placed between the two 
bright central galaxies. 

The Markov chains have a total length of 250 000. After the 
first -10 000 parameter sets w/, the acceptance rate has stabilised 
at 60% as well as the scatter in the x^ values for the parameter 
sets has stabilised. Although a 'burning in' of a Markov chain is 
strictly not necessary, we discard the first 10 000 parameter sets 
as these still show some memory on the initial conditions where 
haloes had no mass. 

Each Markov chain has a fixed description for the mass in 
the galaxies. This means that individual Markov chains do not 
take into account the uncertainties of the velocity dispersion es- 
timates of the galaxies. In order to account also for these uncer- 
tainties we calculate a total of 200 Markov chains where each has 
a diff'erent realisation of the galaxy component. In each Markov 
chain we have reassigned galaxy velocity dispersions from a 
Gaussian distribution centred on the measured value and with a 
width corresponding to the estimated error. Half of the Markov 
chains have a smooth component described by NSIE haloes and 
the other half by ENFW haloes. 

All the remaining analysis is performed by combining the 
200 Markov chains and randomly selecting 1000 parameter sets 
(with replacement) from the combined chain. 



4.0.15. Markov chain Monte Carlo 

The Markov chain Monte Carlo method is a popular way to es- 
timate the errors in model parameters. It allows one to estimate 
the conditional probability distribution of parameters by statis- 
tically sampling the parameter space. In the Metropolis-Hasting 
algorithm the parameter space is explored in a random walk. It is 
a rejection sampling algorithm which means that the next set of 
random parameters is accepted with a probability that depends 
only on the diff'erence in probability between successive sets of 
parameters. Given parameter sets w/ with;^^(w/) and a new per- 
turbed parameter set w/+i with;^^(w/+i) the probability p of ac- 
cepting the new parameter set w/+i is given by the ratio of the 
likelihoods of the parameter sets: 



5. Results 

The multiple images are reproduced well both in position and 
in morphology, although the mean separation between the ob- 
served image position and that predicted by the models is still 
well above the accuracy with which positions of galaxies can 
be measured. The mean separation in the MCMC analysis is 
~ (for reference, the mean separation between an observed 
image position and that predicted by the models is ~ 3.5'' when 
the smooth mass component is described with only one halo). 
For comparison, the best models have separations below ~ 0.5''. 
Observed and predicted morphologies of the images can be seen 
in Figs. lA.ll and lA.2l in the appendix, where a discussion can also 
be found. 
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We explore the mass distribution of the cluster in two dif- 
ferent ways. First of all we can study the parameter space of 
the smooth mass component. Secondly, we look at the total 2- 
dimensional mass distribution of the cluster as well as the aver- 
aged radial mass profile. 



5.1. Parameters of the mass smooth component 

The distribution of the positions of the smooth DM haloes in 
the MCMC analusis are shown in the two top panels in Fig. 
O The left panels are for the NSIE halos and the right pan- 
els for the ENFW haloes. We have used the positions of the 
haloes to separate the primary and secondary haloes. The pri- 
mary halo is defined as the one whose centroid lies within the 
red circle shown in the top panels. In a small number of cases 
both the primary and the secondary halo reside either inside or 
outside the red circle. In these cases we have taken the halo 
with a larger M(r < 34'') as the primary halo. The points are 
coloured according to the mass of the haloes with smaller haloes, 
M(r < 34'0 < 5 X lO^^M©, being light orange and massive 
haloes, M(r < 34'') > 2 x IO^^Mq, being black. It is clear from 
the plot that the central haloes are also the more massive ones, 
and that the less massive haloes populate the outer regions of the 
cluster centre. 

The more massive haloes are spatially concentrated in a very 
small region for both the NSIE and ENFW profiles. The less 
massive ones are more spread out for both profiles, although the 
NSIE profiles are more scattered than those described by ENFW 
profiles. 

In order to be able to compare the masses and concentra- 
tions of the NSIE and ENFW haloes we calculate the mass of 
the haloes at two diff'erent radii ri and r2, where ri < r2. We 
can then define concentration as the ratio of the two masses 

= M(r <r\) I M(r < r2) = Mi / M2. If all the mass is al- 
ready contained within ri then the mass ratio is 1 , where as for a 
mass sheet the ratio is given by (ri/r2)^. For a singular isother- 
mal sphere the maximum value for Cm is (ri/r2), since M(r) oc r. 
The outer radius r2 is taken to be 34 arcseconds which corre- 
sponds roughly to Einstein radius of the cluster, while the inner 
radius is arbitrarily defined as ri = r2/2. The separation in con- 
centration and mass is not very sensitive to the chosen radii, ri 
and r2. On the bottom panels of Fig. [3] we plot the concentration 
Cm against M(r < 34^0 for the NSIE (left) and the ENFW (right) 
haloes. The primary haloes are coloured black, while the sec- 
ondary haloes are coloured grey. For the NSIE profiles the sep- 
aration between the primary and secondary haloes in the Cm vs. 
M(r < 34") space is clear with primary haloes being both more 
massive and more concentrated with a mean mass of M(r < 34") 
= ( 16.1 ± 3.3 ) X 10^^ Mq and a mean concentration of Cm 
= 0.38 ± 0.04 compared to the secondary haloes which have a 
mean mass of M(r < 34'') = ( 8.9 ± 2.8 ) x 10^^ M© and a mean 
concentration of = 0.30 ± 0.06. This separation is less pro- 
nounced for the ENFW profile where the concentrations are vir- 
tually same for the primary and the secondary haloes (c^ = 0.38 
± 0.06 vs. Cm = 0.38 ± 0.11). The mean masses of the primary 
and the secondary haloes are M(r < 34") = ( 16.3 ± 5.0 ) x 10^^ 
Mq and M(r < 34'') = (1.9± 5.4 ) x 10^^ M© respectively. One 
should keep in mind that since the ENFW halo has always an 
inner logarithmic slope of - 1 , the profile can never reach as low 
a concentration as an NSIE profile can. This is also the reason 
why the positions of the secondary haloes for the ENFW profile 
have less scatter: the position of an ENFW halo is always well 
defined. 
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Fig. 5. Mass profile of RX J1347. In the bottom panel we show 
the mass in the galaxies as circles and the total mass as crosses. 
The error bars show the 68% confidence interval derived from 
the MCMC analysis. The top panel shows the residual mass 
needed to perfectly reproduce the image positions as a percent- 
age of the total mass. It is consistent with 0% at all radii. 



5.2. Total mass distribution of RX J1347- 1 145 

Due to the large scatter in the parameters of the smooth mass 
distribution it is perhaps more interesting to look at the total 2- 
dimensional projected mass distribution of the cluster. The av- 
eraged 2-dimensional projected scaled mass distribution for a 
source at redshift 2 is shown in Fig.H] Thick contours lines are 
used for a: < 1 levels, where as thin lines are used for /c < 1 
contours (the contours are separated by Aa:=0.25, the first thick 
contour level is at /c=l). 

For each of the 1000 selected parameter sets from the 
MCMC analysis we have also calculated separately the contribu- 
tions to the total mass from the galaxies and the smooth compo- 
nent. The radial mass profiles of the galaxy mass component and 
that of the total cluster mass are shown in Fig. [5] In the bottom 
panel the total mass (galaxies and the smooth mass component) 
are shown as crosses, and the mass in the galaxies as circles. The 
errors show the 1-sigma errors derived from the MCMC analy- 
sis. The error in the total mass is very small at r ~ 35" where the 
mass is fixed by the spectroscopic redshift of image system 2. 

We can estimate the NSIS and NEW parameters of the total 
mass distribution by fitting a single halo to the total radial mass 
profile. The 1-, 2- and 3-sigma contours for the NSIS and NEW 
parameters are shown in Figs.[6]and|7]respectively. We also plot 
in the figures the error bars for the individual parameters when 
marginalised over the other parameter. The best fit parameters 
of the NSIS profile are o-=1949 k^/s, and rc=20.3 +|;^ ". 
The velocity dispersion is high due to the presence of a large 
core. The NEW profile parameters are r2oo = 3.3 ± 0.2 Mpc, and 
c = 5.4 ^Q^. The;^^ values for both of the profiles are very good 
(for 14 degrees of freedom values of 4.5 and 4.2 are obtained 
for NSIS and NEW profiles respectively). The very good;^^ val- 
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Fig. 3. The distribution of the positions of the two smooth haloes in the MCMC sampling are shown on the two top panels, the 
NSIE profiles are on the left and the ENFW haloes on the right. One of the haloes resides preferentially in a relatively small region 
of the cluster centre, indicated by the red circle. This allows us to make a separation between the primary halo (within the circle) 
and the secondary halo (outside the circle). The points are coloured according to the mass of the haloes. The smaller haloes with 
M(r < 34^0 < 5 X 10 ^^M© are light orange in colour and the massive haloes with M(r < 34'') > 2 x 10 ^^M© are coloured black. In 
the bottom panels we plot the concentration as measured by the ratio of the halo mass at two diff'erent radii, M(r < 17'')/M(r < 34''), 
against the mass within the Einstein radius of the cluster, M(r < 34"). The primary haloes (black) are generally both more massive 
and more concentrated then the corresponding secondary haloes (grey). The separation is clear for the NSIE haloes, where as the 
concentration of the ENFW profiles are the same for both the primary and secondary haloes. The positions on the upper panels are 
relative to = 13:47:30.9, 6 = -11:45:08.6, the centre of the circle. 



ues probably reflect the way the total mass profile was obtained 
using haloes that had NSIE and ENFW profiles. 

5.2.1. RXJ1 347-1 145 as a merger 

It is clear from the surface mass density contours that the mass 
distribution of the cluster does not have any evident bimodality, 
although two haloes were used in modelling the mass. If the clus- 
ter was undergoing a merger along the line of sight, this could 
still require two haloes in the modelling. The fairly clean sepa- 
ration of the two haloes in both mass and concentration as well 
as to a lesser extent also in position, if taken seriously, would 
indicate that RX J 1347 is a merger with a 2:1 mass ratio of 
the merger progenitors. This supports well the results from both 



X-ray s fe.g. lAllen etaDl2002h and SZ-eff'ect (e.g. lKomatsu et al.l 
[2001) that find X-ray shocks and complicated substructure in the 
cluster. This w as also seen in a com bined strong and weak lens- 
ing analysis by lBradac et al.l (l2005l) . Although the 2-dimensional 
mass map obtaned here does not have a clear peak in the south- 
east part of the cluster where the X-ray and SZ structures are 
seen, it is worth noting that it is also in this part of the cluster 
where the second halo is preferentially located. 

We have added a third halo to the smooth mass component 
but have not been able recover the south-east extension. The 
strong lensing constraints in this part of the cluster are weak due 
to the lack of multiple images in this region, and aff'ects our abil- 
ity to determine the mass distribution accurately. This includes 
the position of the secondary halo. 
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Fig. 4. The average scaled surface mass density contours for a source at redshift 2. The thick contours show the k levels greater or 
equal to one, and the thin lines for k less than one (the contours are separated by Aa:=0.25, first thick contour level is at k=1). 



Interestingly the cluster has a very low velocity disper sion of 
only 910+130 km/s measured bv ' Cohen & KneibI (|2002'). They 
discussed that this could be explained if the merger took place 
perpendicular to the line of sight and therefore would have in- 
creased the velocity dispersion mostly in that direction. 

It is important to keep in mind that in the strong lensing anal- 
ysis presented here we assume that the cluster has two haloes. 
More work is necessary in order to establish with more con- 
fidence that the two haloes in our modelling do correspond to 
physical haloes. It is also possible that the need for two haloes in 
the modelling has to do less with physically separated structures 
in the cluster than with a complex morphology of one dominant 
halo. This could be, for example, a varying ellipticity (both mag- 



nitude and position angle) as a function of cluster centric radius, 
slightly asymmetric mass profile, or a mass profile that is neither 
an ENFW nor an NSIE profile but can be modelled as a super- 
position of two such haloes. 

5.2.2. Residual mass maps 

We have taken inspiration from LensPerfect by Dan Coe that, 
as its name implies, allows one to find a mass distribution that 
exactly reproduces all multiple image positions. We have imple- 
mented a similar scheme with the intention to check how our 
parametric mass profile needs to be modified in order to exactly 
reproduce the image positions. The usual parametric modelling 
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Fig. 6. The 1-, 2- and 3-cr confidence regions for a NSIS halo 
fitted to the total mass. The point with error bar shows the best 
fit value and the marginalised 1-sigma errors for the two param- 
eters. The best fit parameters are cr=l949 km/s, and rc=20.3 
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Fig. 7. The 1-, 2- and 3-cr confidence regions for an NFW halo 
fitted to the total mass. The point with error bar shows the best 
fit value and the marginalised 1-sigma errors for the two param- 
eters. The best fit parameters are r2oo = 3.3 ± 0.2 Mpc, and c = 



is used as a basis after which we calculate residual mass map 
needed for a perfect fit. In calculating the residual mass maps 
we have not considered the relative magnifications of the multi- 
ple images but only their positions. 

In practice the residual mass map is found by calculating the 
residual deflection angle Aa(0i j) at each multiple image position 

Aa(6ij) = Pij - (7) 

where Pi j is the source position of image j in system / as ob- 
tained with the parametrised model and is the mean source 
position of the images in system /. 

The residual deflection angle Aa{0) at any other point is then 
calculated using a Thin Plate Spline int erpolation (TPS, for de- 
tails on the topic see e . g . I B ooksteinll 198 9). The interpolated sur- 
face has the property that it minimises the bending energy of the 
function Aa{6) , 

r r\d^Aa I d^Aa f d^Aa 1 , , 
^= — ^+2hr^ + — ^ d6/id^2, (8) 
J J [ de\ \dO1dO2) del \ 

where Oi and 62 are the Cartesian components of the position 
vector 6. 

The amount and location of residual mass Ak(6) that is 
needed to perfectly reproduce the multiple image positions is 
then simply calculated from the residual deflection angle. Since 
the surface mass density is related to the deflection angle via 
K = Va/2, the use of TPS interpolation essentially minimises 
the gradients of Aa:. 

The 2-dimensional distribution of the residual mass Aa: is 
shown in Fig. [H The red and blue contours show positive and 
negative contributions, respectively, while the white line shows 



the Aa: = level. The very small correction needed is repre- 
sented by the very small separation of the contours of only 0.005. 
Although the correction to the mass is small, it is clear from Fig. 
[8] that the correction has a bipolar distribution so that south-east 
and north-west quadrants have generally a mass deficit in the 
parametric modelling while the south-west and north-east quad- 
rants have a slight mass surplus. 

The radial mass profile of the residual mass map is shown 
on the top panel in Fig. [5] as a percentage of the total mass at 
any given radius. In the strong tensing region of the cluster the 
average residual mass at any radius is consistent with 0% with a 
scatter of only around 1%. 

5.3. Lensing redshift of multiple image system 1 

The lensing models allow us to also constrain the redshifts of 
multiple image systems. This can be done by looking at the dis- 
tribution of the redshifts of the multiple image systems in the 
MCMC sampling of the parameter space. The redshift probabil- 
ity density of multiple image system 1 is shown in Fig. [9l The 
dotted line shows the redshift distribution from the lensing and 
the dot-dashed from the photometric redshift analysis, while the 
solid line shows the total distribution. The lensing redshift dis- 
tributions do not show any dependence on the profile used in 
modelling the smooth mass component of the cluster. The lens- 
ing redshift for image system 1 is ziens= 1.90^^3 1. The photo- 
metric redshift was not used as a prior in the lensing modelling, 
and we only use it to compare to the redshift obtained with the 
lensing models. The two estimates for the source redshift agree 
in the redshift range 1.5-2.0 with a combined mean redshift of 
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Fig. 8. The residual surface mass density contours for a source 
at redshift 2. The red and blue contours show positive and neg- 
ative contributions respectively while the white line shows the 
Aa: = level. The contours are separated by 0.005 in k indicat- 
ing that only a very small correction to the mass map is enough 
to perfectly reproduce the multiple image positions. 



5.4. Constraints from all multiple image systems 

All the previous analysis has been based on models that use only 
image systems 1 and 2 to constrain the mass distribution of the 
cluster. We have repeated the analysis using all multiple image 
systems to check how the mass distribution changes in the pres- 
ence of more constraints, and to predict redshifts for the remain- 
ing multiple image systems. The source redshifts are added as 
free parameters to the models without priors from photometric 
redshifts. 

5.4.1. Mass distribution 

The mass distribution obtained with all the multiple images are 
virtually indistinguishable from that obtained with only the two 
first ones. The surface mass density contours for the mass distri- 
bution derived with constraints from all multiple images systems 
are shown in Fig. [TOl The contour levels are the same as those 
used in Fig.H] The differences are most noticeable around some 
of the cluster galaxies. Also the critical lines are only marginally 
affected when all image systems are used. 

The reason for the similarity of the mass maps is most likely 
the free redshifts of the multiple images combined with the ge- 
ometry of the images. With the exception of image systems 9 
and 1 1 all of the image systems are only on one side of the clus- 
ter. The unconstrained redshifts leave the mass scale free while 
the geometries cannot constrain the position of the mass concen- 
tration. All of the information is hence contained in the first two 
image systems which constrain the position of the mass concen- 
tration through the 5 images in image system 1 while the scale 
is fixed by the spectroscopic redshift of image system 2. 
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Fig. 9. The redshift distribution of image system 1 in the MCMC 
sampling of the parameter space. The dotted line shows the 
lensing redshift probability distribution. The photometric red- 
shift probability distrubution is given by the dot-dashed line. 
Combining the two probability densities results in the solid line. 
The lensing redshift probability densities have a clear cutoffs at 
z ~ 1.5 and z ~ 2.5. The lensing redshift for image system 1 is 
^iens=l-90^Q3|, with the best combined estimate being 1.70^q|2. 



5.4.2. Lensing redshifts of multiple image systems 

The redshift probability densities for image systems 1, 6, 8, 9, 
10, 11, and 12 are shown in Fig.[TT]as dotted lines. The redshifts 
of the other multiple image systems are poorly constrained. The 
lower redshift bump on the probability densities for the image 
systems 9 and 10 is a consequence of their modelling. We have 
only required that a critical line should pass between the two im- 
ages in these systems. For a source at low redshift it is the tan- 
gential critical line, and not the radial one, that passes through 
the two images. The dot-dashed line shows the probability distri- 
bution of the photometric redshifts, the solid line the probability 
density combining both the lensing and photometric redshifts 
probability densities. The photometric redshift probability dis- 
tribution of image system 8 has two distinct peaks, one between 
0.6 and 2.0 and another between redshifts 4 and 5. The lensing 
redshift probability distribution for this image system is broad 
but excludes the lower redshift bump seen in the photometric 
redshift probability distribution. The photometric redshift prob- 
ability densities for the other image systems are very broad and 
constrain the redshifts very little. The low redshift bump in im- 
age systems 9 and 10 can be excluded since in those cases the 
image systems would define a tangential critical curve where as 
it is clear from the ACS images that the images are separated by 
a radial critical curve. 

6. Comparison to previous mass estimates 

There are several previous mass estimates for 
RX J1347-1145 in the literature. In this section we com- 
pare the masses from different methods, namely from the 
kinematics of the cluster galaxies. X-ray emission from the 
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Fig. 10. The average scaled surface mass density contours for a 
source at redshift 2. All image systems were used to constrain the 
mass distribution. The thick contours show the k levels greater or 
equal to one, and the thin lines for k less than one (the contours 
are separated by Aa:=0.25, first thick contour level is at k=1), 
these are the same as in Fig.|4]for easier comparison. The diff'er- 
ences are very small, mostly noticeable by small deviations near 
cluster galaxies. 



cluster gas, weak and strong lensing separately (WL & SL, 
analysis respectively) as well as a combined strong and weak 
lensing (SWL). We have converted results in the literature to the 
cosmology used in this work (Hq = 70 km/s/Mpc, Qni=0.3 and 
Qa=0.7), whenever required. For the direct mass comparison 
we have only used literature values where projected masses are 
given. The converted masses are shown in Table[T]and plotted in 
Fig. [El 

The strong lensing mass estimate by lSahu et al.l (1 19981 S98) 
is almost a factor of 4 higher than what we obtain. This is 
probably caused by the low redshift assumed for the arc by 
ISahu et all (1 19981) and their assumption of spherical symmetry. 
The two purely w eak lensing mass estimates have either h igher 
(lKlinget al.l '2QQ5. K05) or lower ( Fischer & TvsonI 1 1997L F97) 
masses than what we obtain with strong lensing. This indicates 
that there is at least no systematic off'set between weak and 
strong lensing masses obtained. 

The agreenient with X-r ay mass estim ates by both 
ISchindler et aP (1 19971 S97) and lGitti et all (^7, G07) is fairly 
good, although all of the X-ray mass estimates are lower than 
the ones we obtain. A detailed comparison of lensing mass from 
a combined strong and weak lensing analysis to new X-ray data 
taken with the Chandra X-ray telescope is presented in Bradac 
et al. (2008, in preparation). We find excellent agreeme nt with 
the co mbined strong and weak lensing analysis of Br adac et al.l 
(I2005L BOS) . This is not surpri sing since the redshift for the arc 
used both in lBradac et al.l ('2005) and this work has not changed, 
thus fixing the mass of the cluster at the position of the arc. 

In cases where we have not found projected mass estimates 
for the cluster we can still compare our results for the parame- 



Fig. 11. The redshift distributions of the image systems for 
which the lensing models are able to constrain the redshift. The 
dotted line shows the lensing redshift probability distribution, 
the dot-dashed line the photometric redshift probability distribu- 
tion and the solid line the combined redshift probability distribu- 
tion. All multiple images are used in the lensing modelling. For 
image systems 9 and 10 the model requirement is that a critical 
curve should pass between the two images. This produces the 
low redshift bump in the redshift probability when the tangen- 
tial critical curve and not the radial critical curve passes between 
the two images. As it is clear from the positions of the multi- 
ple images that they form a radial arc, the lower redshifts can be 
ignored. 

Table 1. Comparison of projected mass estimates interior to a 
given radius in the literature. The diff'erent methods are weak 
lensing (WL), strong lensing (SL), strong and weak lensing 
(SWL), X-rays and kinematics of the cluster galaxies. We have 
converted all the results to the cosmology used in the this paper 
(Ho = 70 km/s/Mpc, Qm =0 -3 and Da =0.7) if necessary. The re- 
marks for lGitti et al.l (|2007|) refer to the types of profiles fitted, 
SO for a single jS model and DDgl for a double j3 model. 
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Gitti et al. ^2007) 35 
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Sahuetal. (1998) 35 

This work 35 

Bradac et al. (2005) 89 

This work 89 

Kling et al. (2005) 124 

Kling et al. (2005) 124 

This work 124 

Gitti et al. (2007) 147 

Gitti et al. (2007) 147 
Fischer & Tvson (1997) 294 
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Fig. 12. Projected mass profile of RX J1347 from various pub- 
lished results. The solid line shows the mass profile obtained in 
this work (the dashed lines show the 1-sigma confidence region 
in mass). The diff'erent symbols in the figure denote diff'erent 
methods, weak lensing (WL), strong lensing (SL), strong and 
weak lensing (SWL) and finally X-rays shown by filled circles, 
crosses, a triangle, and squares respectively (also labelled in the 
figure). The authors are marked next to the points, full details can 
be found in the text. We have plotted only X-ray masses where 
projected masses are given in literature. 



Table 2. Comparison between the estimated parameters of dif- 
ferent parametrised mass profiles from the literature. When com- 
paring parametrised halo profiles there is no confusion between 
spherical and cylindrical densities/masses making the compari- 
son of diff'erent methods easier. We have converted the of an 
isothermal profile and the r2oo of an NFW profile to our cosmol- 
ogy where applicable. If only of an NFW profile was given 
in literature we have multiplied it (and its errors) by the best fit 
concentration to obtain an r2oo- 



Reference 


cr 
km/s 


kpc 


Method 


Cohen & Kneib (2002) 


910 + 130 




kinematic 


Allen et al. (2002) 


1590 + 150 


38 + 8 


X-ray 


Kling eyhQOOS) 
Fischer & Tvson (1997) 


1400 +1^^ 
1500 + 160 




WL 
WL 


This work 


1949 + 40 


117 + 12 


SL 




Reference 


c 


^200 

Mpc 


Method 


Allen et al. (2002) 
Gitti et al. (2007) 


5.87 + 1.4 
3.2 + 0.3 


1 QQ +1.19 
^•^^ -0.60 

2.31 +0.36 


X-ray 
X-ray 


Klins et al. (2005) 
This work 


15 

5 4 +0^^ 

■^'^ -0.5 


2 64 

^•"^ -0.57 

3.29 + 0.20 


WL 
SL 



at the virial radius of the cluster, the value we obtain for r2oo 
is therefore an extrapolation of the NFW profile from the very 
inner regions of the cluster. In our case the scale radius NFW 
profile is already -600 kpc and outside the region that can be 
probed directly with strong lensing only. Weak lensing analysis 
beyond the field of view of the ACS is necessary to strongly 
constrain the r2oo of the cluster. This work will be done in Wuttke 
et al. (2008, in preparation) using both WFI and MegaCam data. 



ters of the commonly used halo profiles, namely the family of 
isothermal profiles and the Navarro, Frenk and White profile. 
These are shown in Table [2l For the comparison we have con- 
verted the literature values of the of an isothermal profile and 
the r2oo of an NFW profile to our cosmology where applicable. In 
addition, if only of an NFW profile was given in the literature 
we have multiplied it (and its error) by the concentration param- 
eter given to obtain an r2oo- This is only a rough estimate of the 
error, however, since the two parameters have very strong degen- 
eracies with higher concentrations having smaller scale radii. 

The high core radius required to fit the mass profile in this 
work is compensated by an unusually high velocity dispersion. 
A singular isothermal sphere would have a much lower velocity 
dispersion (in the range 1450-1600 km/s) as can be seen at least 
qualitatively by extending the degeneracy of the two parameters 
in Fig. [6l The exact value obtained for the velocity dispersion 
depends very strongly on the range of radii in which the mass is 
fitted. On the other hand the range of velocity dispersions is in 
good agreement with the singular isothermal fits in the literature. 
This als o includes the X-ra y based velocity dispersion measure- 
ment bv I Allen et all (l2002h . once the large core radius found in 
our analysis is taken into account (for a given velocity dispersion 
a larger core radius lowers the mass of an isothermal sphere). 

The concentration parameter we obtain for the NFW profile 
is in agreement with the relatively low concentrations expected 
of cluster -sized haloes. The r2oo (3.3 ± 0.2 Mpc) on the other 
hand is rather large when compared to the literature. Since we 
are only determining the mass in the strong lensing regime we 
are not able to eff'ectively constrain the mass at large radii, e.g. 



7. Conclusion 

In this paper we present the first detailed strong lensing model 
for the galaxy cluster RX Jl 347- 1145 . The models are based 
on images taken with the Advanced Camera for Surveys on the 
Hubble Space Telescope. The high resolution of the image along 
with the three colours have allowed us to identify several new 
strongly lensed background images in the cluster. The follow- 
up spectroscopic observations on the F0RS2 at the VLT have 
provided us with redshifts for the images in one multiple image 
system. This is crucial for an accurate absolute mass determina- 
tion of the cluster. Among the new multiple image systems is a 
5 image system with a central image that is very important in 
fixing the mass profile at the core of the cluster. 

The mass in the cluster is modelled with small-scale mass 
components in the cluster galaxies described by truncated 
isothermal ellipsoids, and a large-scale component in two para- 
metric mass profiles (Navarro, Frenk and White profile and non- 
singular isothermal ellipsoid). The parameters of the large-scale 
component are constrained by the location and magnifications 
of the multiple images. The total mass profile of the cluster is 
well described by both a Navarro, Frenk and White profile with 
a moderate concentration of c = 5.3^^^ and r2oo = 3.3^^^ Mpc, or 
a non- singular isothermal sphere with velocity dispersion (Tnsis = 
1949 ± 40 km/s and a core radius - 20.3 ±1.8 " . The core ra- 
dius in the isothermal profile is necessary in order to reproduce 
the central image and to fit the total mass of the profile at the 
all radii. A singular isothermal sphere fit to the total mass has 
too much mass at small radii and too little at larger radii. The 



14 



A. Halkola et al.: The mass distribution of RX J1347-1145 from strong lensing 



total mass of the cluster inside the Einstein radius of the main 
arc (~ 35'^ or -200 kpc) is (2.6 ± 0.12) x 10^^ 

A comparison with the X-ray mass estimates by iGitti et al.l 

(l2QQ7h shows that the lensing mass is consistently higher than 
the X-ray mass, although s till within the erro r bars. The X-ray 
based velocity dispersion in I Allen etal ] (l2002h on the other hand 
is in good agreement with the one determined here once the ef- 
fect of the different core radii is taken into account. 
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Appendix A: Multiple image systems 

The multiple images with some of their properties are listed in 
Table I A. Il The images in each image system are shown in Figs. 
IA.1HA.6[ In addition we show for multiple image systems 1 and 
2 the images and predicted images in Figs. lA.ll and lA.2[ In the 
figures the top row shows the original 3 -colour image created 
from the HST ACS images in F475W (blue), F814W (green) and 
F850LP (red) pass bands. The scales are shown on the images. 
The model prediction on the lower row are created by taking 
the first image of an image system, delensing it into the source 
plane using the lensing model and then relensing it back to the 
image plane at the positions of the multiple images. The first 
multiple image in each system is therefore reproduced exactly 
both in shape and position, while the remaining multiple images 
have an off'set and a distortion relative to the observed multiple 
image due to the model. Note that this assumes that the source 
position of the multiple image system is at the source position 
of the first image, not the average source position which is used 
in calculating the;^^. The sizes and positions are matched well 
indicating that both the magnifications (size) and positions are 
well reproduced by the model. Image lb has somewhat larger 
discrepancy both in magnification and position (the positional 
off'set is 3'0- This can be at least partly due to the presence of a 
cluster galaxy near the multiple image position. A comparison of 
the orientations of the model images with the observed ones (see 
Fig. lA.lb can be used to check the quality of the modelling since 
the orientations have not been used in determining the mass dis- 
tribution. The orentations are in excellent agreement providing 
strong support for the mass modelling presented in this paper. 
The mean separation between an image in an image system and 
that of the models is ~ V in the MCMC analysis. Note that this 
is relatively large since the MCMC chain does not explicitly try 
to minimise the;^^. The best models in the MCMC analysis have 
mean separations below -0.5''. 
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Table A.l. Positions, colours and redshifts estimates for the multiple images in this study. The mass profile and lensing redshifts 
are constrained by image systems 1 and 2 only. The more tentative images in a system are marked with a question mark. 
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J2000 


J2000 










la 


13:47:33.009 
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- 
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44:54.23 


0.12 + 0.03 
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0.806^ 


1.26 + 0.02 


4a 


13:47:31.148 


-11 


44:38.87 


2.59 + 0.03 


-0.76 + 0.01 


0.785^ 


1.32 + 0.02 


5a 


13:47:27.988 


-11 


45:56.48 


1.68 + 0.16 


-1.07 + 0.10 




- 


5b 


13:47:27.720 


-11 


45:51.28 


1.77 + 0.16 


-1.02 + 0.09 




- 


5c 


13:47:27.793 


-11 


45:54.65 


- 


— 




- 


6a 


13:47:29.256 


-11 


45:54.36 


0.55 + 0.12 


-1.63 + 0.25 




0.57 + 0.14 


7a 


13:47:34.797 


-11 


45:01.50 


2.45 + 0.01 


-0.83 + 0.00 




0.44 + 0.06 


8a 


13:47:32.199 


-11 


44:30.48 


0.98 + 0.13 


-1.50 + 0.18 




1.88 + 0.19 


8b 


13:47:31.901 


-11 


44:28.33 


0.83 + 0.10 


-1.73 + 0.18 




- 


8c? 


13:47:34.185 


-11 


45:05.25 


1.00 + 0.11 
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^ Redshift from Lombardi et al. (2008, in preparation). 
2 Redshift from Ravindranath & Ho (2002). 
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Fig. A.l. Image system 1: The offsets between the observed and predicted image positions represent the quality of the models. The 
relatively large offset present in multiple image lb can result from the nearby cluster galaxy. 
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Fig. A.2. Image system 2: 
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Fig. A.3. Image systems 3-7: 
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Fig. A.4. Image systems 8, 9 and 10: 
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Fig. A.5. Image system 1 1 : 
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Fig. A.6. Image systems 12 and 13: 



